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ABSTRACT 

A  static  pressure  probe  was  tested  to  determine  the  feasibility  of 
using  the  probe,  as  an  integral  part  of  a  missile  nose,  to  sense  missile 
altitude.  Experiments  were  conducted  at  Mach  2.0  and  at  Mach  1.51. 

At  Mach  2.0,  the  static  pressure  probe  will  perform  within  altitude 
specifications  of  25,000  feet  ±  2,000  feet  at  angles  of  attack  ranging 
from  -  8  to  +  8  degrees.  At  Mach  2.0,  within  an  angle-of-attack  ranging 
from  0  to  6  degrees,  the  probe  will  measure  free  stream  static  pressure 
within  4  percent;  a  4  percent  error  in  measurement  is  equivalent  to  an 
altitude  error  of  900  feet.  The  missile  nose  shock  will  remain  downstream 
of  the  probe  pressure  ports  for  flight  Mach  numbers  above  1.5. 
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I.   INTRODUCTION 

A.   BACKGROUND 

The  ZEPPO  Rocket  is  an  air-launched  system  designed  to  deliver  an 
expendable  pulse  power  communications  jammer  at  long  range.  To  achieve 
this  objective,  the  system  is  loft-launched  at  angles  up  to  45  degrees 
and  at  high  altitude.  At  the  proper  point  in  the  trajectory,  a  deploy- 
ment sequence  is  initiated  by  the  fuze.  At  the  end  of  this  sequence,  the 
expendable  jammer  is  suspended  and  operated  from  a  17-foot  diameter  para- 
chute with  a  descent  rate  of  10  to  12  feet  per  second.  The  ZEPPO  system 
is  5  inches  in  diameter  and  115  inches  in  length,  and  weighs  approximately 
135  pounds.  The  deployment  system  consists  of:  (1)  a  timer  fuze  in  the 
nose;  (2)  an  expulsion  charge;  (3)  a  drogue  parachute;  (4)  a  main  para- 
chute; and  (5)  an  expendable  electronic  jammer. 

The  rocket  motor  is  the  production  MK71  MOD  1  Zuni  motor.  The  motor 
weighs  80  pounds,  and  is  73.3  inches  in  length  and  5  inches  in  diameter. 
It  has  a  thrust  level  of  6500  pounds  force  for  1.5  seconds,  with  a  burn- 
out spin  rate  of  30  to  33  revolutions  per  second. 

The  fuze,  at  the  present  time,  is  a  modified  MK375  MOD  0  timer  fuze 
with  a  timing  capability  up  to  140  seconds.  It  contains  a  black  powder  ex- 
pulsion charge. 

The  drogue  parachute  has  a  cruciform  geometry  and  is  packed  in  a 
small  bag  that  opens  after  motor  separation.  The  main  parachute  is  a 
lightweight  17-foot  parachute  adapted  from  the  LUU-2/B  flare  system  which 
is  a  standard  production  item.  The  chute  is  packed  in  a  deployment  sleeve 
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measuring  5  inches  in  diameter  by  eight  feet  in  length;  this  sleeve  dras- 
tically reduces  snatch  and  opening  loads,  permitting  the  lightweight  para- 
chute to  survive  high  "Q"  deployment  without  damage. 

The  optimum  point  for  deployment  initiation  is  25,000  feet,  plus  or 
minus  2,000  feet.  Deployment  at  this  altitude  allows  maximum  payload  time 
aloft,  and  therefore  maximum  operation  time.  System  testing  by  the  Naval 
Weapons  Center,  China  Lake,  California,  has  demonstrated  that  there  is 
difficulty  in  attaining  the  optimum  deployment  altitude  with  the  present 
fuzing  system. 

To  remedy  this  situation,  several  alternatives  were  considered,  one 
of  which  was  a  fuzing  system  initiated  by  a  pressure-sensing  switch,  or 
Baro-Fuze.  Two  alternatives  were  considered  for  measuring  pressure:  a 
shock-swallowing  air  sensor  and  a  static-pressure  probe.  The  Naval  Post- 
graduate School  was  selected  to  evaluate  the  effectiveness  of  the  static- 
pressure  probe.  The  probe  selected  for  this  evaluation  was  developed  by 
Pinckney  [1] . 

Figure  1-1  shows  the  probe  interior  to  the  missile  nose,  where  it 
remains  until  after  launch.  After  rocket  motor  burnout,  and  with  the  mis- 
sile flying  at  small  angles -of-attack,  the  probe  is  extended.  With  the 
probe  extended,  as  shown  in  Figure  1-2,  the  aerodynamic  drag  is  reduced; 
thus  the  decay  of  missile  velocity  is  decreased.  The  range  of  possible 
Mach  numbers  was  considered  to  be  Mach  1.4  to  Mach  2.0  at  small  angles - 
of-attack. 
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B.   THESIS  PURPOSE 

The  purposes  of  this  thesis  are:  (1)  to  evaluate  theoretically  the 
static  pressure  probe  general  designs  developed  by  Pinckney  [1]  for  use 
in  the  Mach  1.4  to  Mach  2.0  range;  (2)  to  develop  pressure  distributions 
along  the  probe,  and  select  the  probe  design  best  suited  to  this  Mach 
number  range;  (3)  to  evaluate  experimentally  the  probe  pressure-sensing 
performance  as  an  integral  part  of  the  missile  nose  (with  probe  extended) ; 
and  (4)  to  report  on  experimentation  in  the  Naval  Postgraduate  School 
supersonic  wind  tunnel  utilizing  the  existing  Mach  1.4  and  Mach  2.0  noz- 
zle blocks. 
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II.   THEORY 

A.   STATIC  PRESSURE  PROBE  DESIGN 

Pinckney  [1]  discusses  the  design,  for  supersonic  flow,  of  a  short 
static  pressure  probe  that  is  relatively  insensitive  to  angle-of -attack. 
The  general  probe  design,  which  is  illustrated  in  Figure  2-1,  consists  of 
a  conical  nose  followed  by  a  parabolic  tangent  ogive  fairing  the  nose  into 
a  second,  truncated  cone.  The  second  cone  mates  with  the  cylindrical  after- 
body of  the  probe.  The  parabolic  formula  used  for  the  tangent  ogive  sec- 
tion was  of  the  form: 


I  =  ±  [A(|  )2  +  B(|  )  +  C]h  (2-1) 

n  n  n 


The  value  of  Y/R  at  the  junction  between  the  tangent  ogive  and  the 
second  truncated  cone  was  specified.  The  conical  nose  half -angle  (3)  and 
the  second  truncated  cone  half-angle  (u)  were  also  specified.  Based  on 
this  information,  the  constants  A,  B,  and  C  were  calculated.  The  probe 
nose  half -angle  is  dependent  on  the  free-stream  Mach  number,  M^;  the  cone 
angle  must  be  small  enough  so  that  shock  detachment  does  not  occur.  Static 
pressure  holes  are  located  at  points  where  static  pressure  equals  free- 
stream  static  pressure.  These  points  of  pressure  equality  vary  with  the 
length  and  angle  of  the  conical  nose.  Pinckney  [1]  calculated  the  theor- 
etical pressure  distribution  for  this  general  design  utilizing  a  blunt- 
body  program  coupled  with  a  method  of  characteristics  program.  He  made 
calculations  for  3  =  10  degrees  with  u  varying  between  2  and  4  degrees 
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for  MOT  =  2.5  and  4.0  and  zero  angle-of-attack.  Wind  tunnel  tests  were  con- 
ducted utilizing  four  different  probe  types.  The  first  two  probe  types  had 
S  =  10  degrees  and  oj  =  2  or  3  degrees,  and  were  tested  at  M,,,  =  2.5  and  4. 
The  second  two  probe  types,  which  had  6  =  20  degrees  and  oj  =  2.0  or  3.5 
degrees,  were  also  tested  at  H,,  =  2.5  and  4.  Tests  were  conducted  at  attack 
angles  ranging  from  0  to  12  degrees. 

Figure  2-2,  which  is  reproduced  from  Kuethe  and  Chow  [2],  is  a  graph 
of  shock  angle  versus  Mach  number  for  various  cone  angles.  Since  the  thrust 
of  this  research  concerns  Mach  numbers  between  1.4  and  2.0,  Figure  2-2 
shows  that  8  cannot  be  greater  than  approximately  15  degrees,  or  shock  wave 
detachment  will  result  at  MOT  =  1.4.  Therefore,  a  conical  nose  half-angle 
of  6  =  10  degrees  was  chosen.  Pinckney's  theoretical  results  for  8  =  10 
degrees,  which  are  shown  in  Figure  2-3,  indicate  that  the  pressure  distri- 
bution which  most  closely  approximates  that  of  the  free-stream  static  pres- 
sure will  be  obtained  with  oj  =  3  degrees  at  M»  =  2.5.  Pinckney's  experimental 
results,  which  are  shown  in  Figure  2-4,  indicate  close  agreement  with  the 
theoretical  calculations.  Figure  2-5  illustrates  Pinckney's  probe  design 
for  8  =  10  degrees  and  u  =  3  degrees.  Since  Pinckney's  theoretical  calcu- 
lations were  for  much  higher  Mach  numbers,  it  was  considered  prudent,  for 
research  completeness,  to  obtain  theoretical  values  for  the  static  pressure 
distribution  at  lower  Mach  numbers.  For  consistency,  attempts  were  made  to 
obtain  the  computer  program  used  by  Pinckney  [1] ;  however,  the  program  was 
not  immediately  available.  Therefore,  an  alternative  approach  was  sought. 
The  coefficients  of  the  parabolic  tangent  ogive  sections  for  8  =  10  degrees 
and  oj  =  3  degrees  were  readily  available  and  were  provided  by  Mr.  Pinckney. 
The  values  of  the  coefficients  are  as  follows:  A  =  -0.009267312,  B  =  8.3676155, 

16 


and  C  =  -402.88107.  The  coefficients  were  verified,  and  Eq.  2-1  provides 
the  ogive  section  body  shape  used  in  the  theoretical  work  described  in  the 
next  section. 

B.  STATIC  PRESSURE  DISTRIBUTION 

Gawain  and  Schonberger  [3]  discuss  small  perturbation  supersonic  flow 
theory  applied  to  a  body  of  revolution.  The  theory  presented  in  Ref .  3  is 
essentially  that  of  classical  slender-body  small -perturbation  theory,  with 
corrected  boundary  conditions  in  the  sense  that  the  boundary  conditions 
are  applied  at  the  body  surface.  In  more  elementary  theory,  the  boundary 
conditions  are  applied  along  the  axis  of  the  body.  The  use  of  the  exact 
boundary  conditions  yields  more  accurate  solutions  than  those  obtained 
from  more  elementary  theory.  The  theory  of  Ref.  3  was  considered  adequate 
for  the  purposes  of  this  thesis,  as  only  small  angles-of-attack  were  con- 
sidered. A  BASIC  language  computer  program,  suitable  for  use  on  a  micro- 
processor, was  developed  based  on  the  theory  and  calculation  procedures 
in  Ref.  3.  The  program,  and  a  sample  data  run  for  MOT  =  2.0  at  angle-of- 
attack  a  =  0  and  10  degrees,  are  given  in  Appendix  A. 

C.  RESULTS  OF  THEORETICAL  CALCULATIONS 

1.  Theoretical  Calculations 

The  pressure  distribution  was  calculated  for  M^  =  1.4,  1.6,  1.8, 
2.0,  and  2.5  at  zero  angle-of-attack  (a).  The  M^  =  2.5  calculations  served 
as  a  comparison  between  Pinckney's  theoretical  and  experimental  results  [1] 
and  the  theory  used  in  this  research.  Figure  2-6  is  a  graph  of  P/P,*,  versus 
X/P^  for  Mo,  =  2.5,  3  =  10  degrees,  and  u  =  3  degrees.  Comparing  Figures 
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2-3  and  2-6,  it  can  be  seen  that  the  graph  shapes  are  the  same;  however, 
the  theoretical  results  depicted  in  Figure  2-6  are  slightly  higher  than 
those  of  Figure  2-3  (approximately  2.7  percent).  Yet,  in  comparing  Figure 
2-6  with  the  experimental  results  of  Figure  2-4,  only  a  1  percent  differ- 
ence exists. 

In  addition  to  the  calculations  above,  computations  were  made  for 
angles-of-attack  a=l,  2,  4,  6,  8,  and  10  degrees  at  M^  =  1.4,  1.6,  1.8, 
and  2.0.  The  calculation  procedure  and  results  will  be  discussed  in  Sec- 
tion 4,  below. 

2.  Probe  Pressure  Hole  Location 

As  was  discussed  previously,  in  Part  A  of  this  chapter,  optimum 
probe  pressure  hole  location  depends  on  the  free-stream  Mach  number.  Due 
to  the  range  of  Mach  numbers,  a  design  compromise  is  made  that  minimizes 
the  pressure  error  over  the  range  of  M^,  rather  than  for  one  specific  M,,. 
Figure  2-7  is  a  graph  of  P/P^  versus  X/P^  for  M,,,  =  1.4,  1.6,  1.8,  and  2.0, 
at  zero  angle-of-attack.  The  graph  shows  clearly  that  the  minimum  differ- 
ence in  P/Pa,  for  the  range  of  Mo,  occurs  at  approximately  X/R^  =  260.  The 
difference  in  P/POT  due  to  variable  M^  is  approximately  4  percent. 

3.  Error  in  Altitude  Due  to  Error  in  Pressure  Measurement 
Pressure  as  a  function  of  altitude  is  given  by: 

P  =  PQe<h        ,  (2-2) 

where  E,   =  (7.6  km)   and  P  is  the  pressure  at  sea  level.  Eq.  2-2  was  ob- 
tained by  graphing  pressure  versus  altitude  using  data  from  the  U.   S. 
Standard  Atmosphere,    1976   [4],  and  fitting  a  curve  to  the  data.  Taking  the 

18 


natural  logarithm  of  both  sides  of  Eq.  2-2  yields 


In?   =  -01   +  In?  .  (2-3) 


Since  P  =  1  atmosphere,  Eq.  2-3  becomes 


In?   =  -£fe    .  (2-4) 

To  find  the  error  in  altitude  due  to  the  difference  in  pressure,  the  deri- 
vative of  Eq.  2-4  is  taken,  which  gives: 

dV/V   =  -Idh         .  (2-5) 

Rearranging  Eq.  2-5,  the  altitude  error  becomes: 


dh  -   -  i  •  d?/?  .  (2-6) 


With  1/5  =  7.62  km  and  d?/?   =  0.04,  the  value  of  dh   is  300  meters,  or  984 
feet.  This  error  is  well  within  the  tolerance  specified  for  probe  perform- 
ance, i.e.,  25,000  feet  ±  2,000  feet. 

4.  Probe  Pressure  Readings  at  Angle-of-Attack 

Figure  2-8  shows  a  cross  section  of  the  static  pressure  probe. 
The  static  pressure  hole  on  the  windward  side  is  number  1,  while  the  hole 
on  the  leeward  side  is  number  3.  Using  the  principle  of  conservation  of 
mass,  the  mass  flow  rate  in,  m . ,  must  equal  the  mass  flow  rate  out,  m   ; 
therefore: 
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PjVjA,  ♦  2p2V2A2  *  p3V3A3  =  0    .  (2-7) 


Since  p.  =  p  =  p  and  A.  =  A2  =  A, ,  Eq.  2-7  reduces  to 


Vj  +  2V2  +  V3  =  0    .  (2-8) 


Note  that  Eqs.  2-7  and  2-8  require  the  velocity  to  be  an  algebraic  quan- 
tity: a  flow  into  the  probe  is  defined  as  positive.  Streeter  and  Wylie  [5] 
demonstrate  that  flow  through  an  orifice  can  be  represented  by: 


PQ  -  P  =  %CpV2    ,  (2-9) 

where  PQ  is  the  pressure  upstream  of  the  orifice,  P  is  the  pressure  down- 
stream of  the  orifice,  and  C  is  a  constant  for  incompressible  flow.  Writing 
Eq.  2-9  in  terms  of  the  parameters  of  Figure  2-8  yields: 

lPl  '  Ppl  =  ^piV?Ci  (2"10) 

lP2"Ppl  =^2V2C2  (2-U) 

lP3  -?pl  =^3V?C3     •  (2"12} 

where  P  is  the  pressure  measured  by  the  probe  and  C  is  the  flow  coeffici- 
ent for  pressure  drop  across  the  pressure  holes.  Solving  Eqs.  2-10,  2-11, 
and  2-12  for  Vj ,  V2,  and  V3,  respectively,  and  substituting  into  Eq.  2-8, 
yields: 

[2  f  P  -P    Ip  ±   2  [2  I  P  -P    \p  ±    [2  |  P   -P    11^=0  .  (2-13) 

L     '     1      p  LI2p|J  L|3p' 

Simplification  to  this  form  was  made  possible  because  p  =  p  =  p  and 
C  =  C  =  C  .  Note  that  the  value  of  C  is  not  needed,  since  C  cancels  in 
Eq.  2-13. 
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Eq.  2-13  is  nonlinear,  and  its  eight  mathematically  possible  cases 
are  shown  in  Table  2.1.  A  plus  sign  (+)  indicates  flow  into  the  probe;  a 
minus  sign  (-)  indicates  flow  out  of  the  probe.  Some  cases  which  are  mathe- 
matically possible  but  not  physically  possible  can  be  eliminated. 


Case  PrPp  P2-Pp     P3-Pp 

1  +  +  + 

2  -  +  + 

3  +  -  + 

4  +  + 

5  + 

6  -  + 

7  -  -  + 
8 

Table  2.1     Mathematically  Possible  Cases  for  Eq.    2-13 


Cases  1  and  8  can  be  eliminated  immediately  as  not  satisfying  con- 
servation of  mass.  In  addition,  at  positive  angles-of-attack  (i.e.,  probe 
tip  up),  the  pressure  P.  will  be  higher  than  Po  or  P  .  Therefore,  cases  2, 
6,  and  7  may  be  ruled  out  as  not  possible.  An  iterative  computer  program 
was  written  to  solve  Eq.  2-13  for  P  ,  given  that  cases  3,  4,  and  5  are 
the  only  possible  cases.  A  listing  of  the  program  and  a  sample  data  run 
can  be  found  in  Appendix  B. 

Figure  2-9,  which  plots  P/P^  as  a  function  of  a,  shows  probe  theo- 
retical performance  for  values  of  a  ranging  between  zero  and  ten  degrees. 
The  data  in  Figure  2-9  were  calculated  for  MOT  =  1.4  and  for  MOT  =  2.0. 
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III.   EXPERIMENTAL  APPARATUS 

A.   WIND  TUNNEL 

1.  Mach  2.0  Wind  Tunnel  Tests 

Wind  tunnel  tests  were  conducted  in  the  Naval  Postgraduate  School 
4x4  inch  supersonic  wind  tunnel  using  the  existing  Mach  2.0  nozzle  blocks. 
(The  Mach  number  of  2  is  nominal;  the  actual  Mach  number  is  1.92).  Figures 
3-1  and  3-2  are  two  photographs  of  the  general  experimental  apparatus.  They 
present  two  views  of  the  supersonic  wind  tunnel  and  its  associated  instru- 
mentation. Utilizing  a  one-fifth  scale  model,  and  regulating  wind  tunnel 
stagnation  pressure,  an  actual  Reynolds  number  of  Re  =  1.3  x  105  was  at- 
tainable, which  compares  favorably  with  the  Reynolds  number  for  the  full- 
scale  probe  at  25,000  feet  (viz.,  Re  =  4.3  x  105) .  The  reference  length 
for  the  calculation  of  the  Reynolds  number  is  0.020  meters.  Twenty-three 
wind  tunnel  tests  were  conducted  at  H*,  =  1.92  with  a  =  0,  1,  2,  4,  6,  8, 
and  10  degrees. 

Prior  to  any  testing,  the  fine-mesh  wire  screen  in  the  wind  tun- 
nel plenum  was  discovered  to  be  torn,  and  hanging  loose  in  the  plenum 
chamber.  The  screen  serves  to  insure  a  flat  velocity  profile  (flat  stag- 
nation pressure  profile)  and  to  generate  small-scale  turbulence  which  dis- 
sipates prior  to  arrival  at  the  test  section.  The  plenum  chamber  exterior 
can  be  seen  in  Figure  3-1.  Sufficient  manpower  to  replace  the  plenum  cham- 
ber screen  was  not  available.  Therefore,  the  screen  was  removed.  To  verify 
that  the  flow  in  the  wind  tunnel  test  section  was  uniform,  an  additional 
pressure  tap  was  installed  opposite  one  of  the  existing  pressure  taps. 
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Testing  showed  no  significant  pressure  difference  (approximately  0.5  per- 
cent) between  the  two  pressure  taps.  Therefore,  the  pressure  measurements 
suggest  that  the  test  section  velocity  profile  was  reasonably  uniform  at 
Mach  2.0. 

2.  Mach  1.4  Wind  Tunnel  Tests 

Tests  were  attempted  utilizing  the  existing  Mach  1.4  nozzle  blocks. 
Supersonic  flow  could  not  be  established  in  the  test  section.  Appendix  C 
provides  a  sequential  listing  of  testing  attempts  and  corrective  measures 
taken  between  tests,  and  discusses  the  three  most  probable  reasons  for 
testing  failure. 

3.  Mach  1.51  Wind  Tunnel  Tests 

Since  experimental  data  at  a  lower  Mach  number  were  still  desired, 
various  alternatives  were  explored.  The  four  alternatives  investigated  were: 
(1)  to  replace  the  plenum  chamber  screen;  (2)  to  design  and  manufacture 
Mach  1.5  nozzle  blocks;  (3)  to  enlarge  the  existing  Mach  1.4  nozzle  block 
second  throat  by  machining;  or  (4)  to  attempts  tests  using  a  combination 
of  Mach  1.4  and  Mach  2.0  nozzle  blocks.  Alternatives  1  through  3  were  ruled 
out  due  to  inadequate  manpower,  funding,  and/or  research  time.  Alternative 
4  was  attempted.  Testing  in  this  unorthodox  manner  was  attempted  in  order 
to  verify  that  the  missile  nose  shock  did  not  detach  at  a  lower  Mach  num- 
ber. Two  shock  waves  occur  at  the  nose  of  the  missile  with  the  probe  extended. 
A  nearly  conical  shock  is  attached  to  the  probe.  A  curved  shock  wave  origi- 
nates due  to  the  missile  ogive;  this  curved  shock  wave  is  termed  the  missile 
nose  shock.  The  probe  conical  shock  and  missile  nose  shock  merge  to  form 
the  bow  shock  wave.  The  probe  conical  shock,  missile  nose  shock,  and  re- 
sulting bow  shock  wave  can  be  seen  in  Figure  3-3,  which  is  a  photograph  of 
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the  model  mounted  in  the  wind  tunnel.  The  probe  is  oriented  at  zero  angle- 
of-attack,  and  the  flow  is  at  MOT  =  1.92.  Both  the  probe  conical  shock  wave 
and  the  missile  nose  shock  wave  are  attached.  If  the  missile  nose  shock 
moves  upstream  of  the  probe  pressure  ports,  erroneous  pressure  readings 
are  obtained. 

The  Mach  number  gradient  in  the  wind  tunnel  test  section  resulting 
from  the  combination  of  nozzle  blocks  was  calculated  using  the  method  of 
characteristics  found  in  Liepmann  and  Roshko  [6] .  Figure  3-4  shows  the 
Mach  number  gradient  in  the  test  section  and  the  shock  waves  that  result 
using  the  combination  of  nozzle  blocks.  Figure  3-5  is  a  photograph  of  the 
Mach  2.0  and  Mach  1.4  nozzle  blocks  installed  in  the  wind  tunnel.  The  theo- 
retical Mach  number  obtained  using  this  combination  of  nozzle  blocks  was 
M^  =  1.51  at  the  probe  nose,  M^  =  1.538  at  the  top,  and  M^  =  1.16  at  the 
bottom  of  the  wind  tunnel  test  section.  Experimental  results,  which  are 
discussed  in  Chapter  IV,  were  used  to  verify  theory. 

B.   MODEL  DESIGN 

A  one-fifth  scale  model  of  the  static  pressure  probe  and  missile  nose 
was  used.  Since  a  complete  model  could  not  be  mounted  in  the  tunnel,  the 
missile  nose  was  truncated  at  40  percent  of  the  nose  length.  Maximum  model 
size  was  dictated  by  maximum  allowed  wind  tunnel  blockage  (approximately 
6.4  percent)  at  Mach  1.4.  The  allowed  wind  tunnel  blockage  limited  the 
combined  frontal  area  of  the  model  and  mounting  system.  Pinckney  [1]  sup- 
plied the  drawings  for  the  static  pressure  probe.  The  Naval  Weapons  Center, 
China  Lake,  California,  provided  the  formula  and  coefficients  for  the  Von 
Karmen  ogive  missile  nose.  The  Von  Karmen  ogive  is  given  by: 
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Y  =  4=  •  [coa"l(l  -  ~5  -  iisinC2cos'1(l  -  ^-))]h  (3-1) 


for  the  full-scale  missile,  where  R  =  2.5  inches  and  L  =  10  inches.  Figure 
3-6  is  a  front  view  of  the  model  and  double-wedge  airfoil  used  to  mount 
the  model  in  the  wind  tunnel.  Design  compromises  between  the  maximum  mis- 
sile nose  diameter  and  maximum  wedge  thickness  were  required  to  remain 
within  wind  tunnel  blockage  limits.  The  maximum  diameter  of  the  model  nose 
used  is  0.792  inches,  and  the  maximum  wedge  thickness  is  0.165  inches.  The 
mounting  wedge  length  is  3.900  inches,  and  a  wedge  half-angle  of  9  degrees 
was  utilized  to  prevent  oblique  shock  detachment  at  Mach  1.4.  The  frontal 
area  of  the  model  plus  mounting  system  was  1.005  square  inches,  which  is 
within  the  allowed  wind  tunnel  blockage.  Figure  3-7  is  a  photograph  of 
the  model  ogive  nose  and  the  extended  static  pressure  probe.  The  ogive  is 
1.200  inches  in  length.  The  pressure  probe  extends  0.680  inches  from  the 
missile  nose.  Figure  3-8  is  a  photograph  of  the  static  pressure  probe 
mated  to  the  missile  nose.  The  four  probe  pressure  holes  are  90  degrees 
apart  and  0.0135  inches  in  diameter,  and  are  located  0.1325  inches  from 
the  probe  conical  nose  tip.  The  parabolic  tangent  ogive  and  second  trun- 
cated cone  start  at  0.0525  inch  and  0.1225  inch,  respectively,  from  the 
probe  nose  tip,  and  are  made  of  0.060-inch  outside  diameter,  0.040-inch 
inside  diameter  stainless  steel  tubing.  The  probe  nose  cone  through  the 
second  truncated  cone  is  made  of  0.0625-inch  diameter  hole  on  centerline. 
After  machining,  the  drill  rod  is  mated  to  the  stainless  steel  tubing 
with  silver  solder.  Figure  3-9  is  a  photograph  of  the  top  view  of  the 
model  and  mounting  wedge.  Figure  3-10  is  a  photograph  of  the  model  mounted 
in  the  wind  tunnel,  and  Figure  3-11  is  a  close-up  of  the  model  mounted  in 

the  wind  tunnel  test  section. 
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C .   INSTRUMENTATION 

Pressure  taps  in  the  wind  tunnel  wall  and  on  the  model  surface  were 
connected  to  a  Giannini  12-position  scanivalve  with  Tygon  tubing.  The 
scanivalve  has  a  range  from  absolute  zero  to  100  PSIA.  The  scanivalve  was 
connected  to  a  Gould  Statham  absolute  pressure  transducer,  Model  PA732TC- 
100-350.  The  transducer  was  powered  by  a  five-volt  excitation;  it  has  a 
pressure  range  from  absolute  zero  to  100  PSIA,  with  a  maximum  frequency 
response  of  8,700  Hz.  A  variable-gain  differential  amplifier  conditioned 
the  transducer  output  to  an  Analog  Devices  4^-digit  voltmeter.  The  voltage 
displayed  is  directly  proportional  to  the  pressure.  Transducer  calibration 
was  accomplished  utilizing  a  60-inch  mercury  column  and  vacuum  pump.  For 
higher  pressure,  a  Wallace  and  Tiernan  bourdon  tube  absolute  pressure  gauge 
was  used  (range  absolute  zero  to  100  PSIA) .  Schlieren  photographs  of  the 
flow,  including  shock  waves,  were  taken.  A  1,000-watt  mercury  vapor  light 
source  and  a  reflex-type  camera  converted  to  accept  Polaroid  Type  152  film 
were  used. 
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IV.   EXPERIMENTAL  RESULTS 

A.   MACH  2.0  EXPERIMENTAL  RESULTS 
1.  Data  Reduction  Technique 

Twenty- three  wind  tunnel  tests  were  conducted  at  Mach  2.0  with 
angles-of-attack  varying  from  -  2  to  +  10  degrees.  Tables  4.1  and  4.2 
depict  sample  wind  tunnel  test  data  for  a   =  0  and  2  degrees,  respectively 

Run:  4    P  =  29.85  in  Hg    Date:  24  September  1980 


=   1.92 

a   =   0° 

T  =   16.2°C 

Time:   0902 

Model 

Port  U 

Plenum 

Reading 

CPSIG) 

CPSIG) 

(PSIG) 

1 

-9.37 

-9.25 

+23.12 

2 

-9.37 

-9.50 

+22.82 

3 

-9.63 

-9.49 

+22.14 

4 

-9.77 

-9.77 

+19.92 

5 

-9.86 

-9.82 

+19.27 

Table  4.1     Sample  Wind  Tunnel  Data  for   a  =  0  Degrees  and   M^  =  1.92 


Run:  6    P  =  29.88  in  Hg    Date:  24  September  1980 
M„  =  1.92 


Reading 

1 
2 
3 

4 
5 

Table  4.2     Sample  Wind  Tunnel  Data  for   a  =  2  Degrees  and   M*,  =  1.92 
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=  2° 

T  =   20.5°C 

Time:    1027 

Model 

Port  §1 

Plenum 

(PSIG) 

(PSIG) 

(PSIG) 

-9.12 

-9.02 

+24.75 

-9.16 

-9.06 

+24.58 

-9.18 

-9.10 

+24.57 

-9.21 

-9.18 

+24.14 

-9.55 

-9.28 

+25.90 

The  ratio  of  measured  static  probe  pressure  to  wind  tunnel  wall 
pressure,  P/POT,  was  calculated  for  each  set  of  data  points.  The  arithmetic 
mean,  variance  (N  weighting),  and  standard  deviation  (N  weighting)  were 
calculated  using  the  calculation  procedure  found  in  Haber  and  Runyon  [7] . 
Table  4.3  gives  these  values  for  various  angles-of-attack,  a. 


Angle-of- 
Attaak 

Theoretical 

P/Pco 

Experimental 
P/P*, 

Standar 
Deviatic 

0° 

1.003 

-  0.987 

±0.00861 

1° 

1.0015 

0.984 

±0.0102 

2° 

0.996 

0.987 

±0.0097 

4° 

0.976 

0.978 

±0.0122 

6° 

0.952 

0.974 

±0.008 

8° 

0.938 

0.932 

±0.013 

10° 

0.918 

0.927 

±0.016 

Table  4.3     Angles-of-Attack  Versus  Theoretical  and 
Experimental   P/P^  for   M^  =  2.0 


2.   Static  Pressure  Probe  Performance 

Figure  4-1  is  a  plot  of  the  data  listed  in  Table  4.3.  Calculated 
points  are  indicated  by  the  letter  X,  while  experimental  mean  values  are 
denoted  by  an  asterisk  (*)  and  standard  deviations  are  denoted  by  the 
character  "#".  The  scale  of  the  ordinate  emphasizes  greatly  any  deviation 
from  P/Poo  =  1.  Note  that  an  experimental  error  bar  is  not  shown  for  angle- 
of-attack;  such  an  error  bar  would  be  horizontal.  Pressure  measurements 
for  angles-of-attack  of  -  1  and  -  2  degrees  are  included  in  the  measure- 
ments shown  for  a  =  1  and  2  degrees .  A  change  in  a  from  +  2  to  -  2  degrees 
had  little  or  no  effect  on  pressure  measurements.  Figure  4-3  shows  that 
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probe  performance  is  reasonably  constant  in  the  range  a  =  0  to  6  degrees; 
the  corresponding  values  of  mean  P/POT  are  0.987  and  0.974.  For  angles-of- 
attack  greater  than  6  degrees,  the  measured  pressure  ratio  is  1.6  percent 
less  than  that  predicted  by  the  theory.  The  mean  value  of  the  measured 
pressure  ratio  is  almost  exactly  that  of  the  theory  at  a  =  4  degrees;  at 
a  =  6  degrees,  the  measured  P/P^  is  2.6  percent  greater  than  theoretical 
predictions. 

Using  Eq.  2-6  with  dh   =  2,000  feet  (609.6  meters)  and  1/?  =  7.62 
kilometers,  the  fractional  change  in  pressure  is  0.08.  Therefore,  taking 
experimental  uncertainty  into  consideration,  the  probe  will  perform  within 
altitude  specifications  for  angles-of-attack  less  than  or  equal  to  8  de- 
grees at  Mach  2.  As  stated  in  Chapter  III,  Figure  3-3  is  a  photograph  of 
the  model  mounted  in  the  wind  tunnel.  The  probe  is  oriented  at  zero  angle- 
of-attack,  and  the  flow  is  at  MOT  =  1.92.  Both  the  probe  conical  shock  wave 
and  the  missile  shock  wave  are  visible  and  attached.  The  horizontal  lines 
just  above  the  model  are  due  to  condensation  within  the  boundary  layer  on 
the  inside  of  the  test  section  window.  The  dark  vertical  line  aft  of  and 
below  the  model  is  the  shadow  of  the  Tygon  tubing  used  to  measure  probe 
pressure.  Figure  4-2  is  a  photograph  of  the  model  oriented  at  a  =  +  1 
degree  within  the  wind  tunnel  operating  at  M^  =  1.92. 

B.   MACH  1.51  EXPERIMENTAL  RESULTS 

Seven  wind  tunnel  tests  were  conducted  using  the  combination  of  Mach 
1.4  and  Mach  2.0  nozzle  blocks.  The  model  and  plenum  pressures,  as  well  as 
pressures  at  five  static  pressure  ports  located  along  the  wind  tunnel  wall, 
were  recorded  for  each  wind  tunnel  test.  The  static  pressure  ports  extended 
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from  just  aft  of  the  nozzle  throat  to  the  test  section.  The  position  of 
the  static  pressure  ports  relative  to  the  wind  tunnel  test  section  is 
shown  in  Figure  3-1.  The  static  pressure  measurements  were  used  to  verify 
the  calculations  based  on  the  method  of  characteristics  referred  to  in 
Chapter  III.  The  measured  static  pressure  for  the  seven  wind  tunnel  tests 
was  recorded,  grouped  by  test  point,  and  the  arithmetic  mean  and  standard 
deviation  were  calculated  for  each  point.  Theoretical  values  and  actual 
experimental  values  were  compared,  and  are  shown  in  Table  4.4. 


Theoretical 

Observed  Mean 

Test 

Pressure 

Pressure 

Theoretical 

Point 

(PSIA) 

(PSIA) 

Mach  Number 

P 

1 

5.97 

5.60 

1.57 

P 

2 

6.05 

6.35 

1.56 

P3 

7.41 

7.92 

1.42 

P 

8.87 

9.12 

1.29 

P 

10.00 

10.89 

1.20 

Percent 
Pressure 

Difference 

5.9 
4.7 
6.5 
2.7 
8.1 


Table  4.4     Theoretical  and  Observed  Pressure  for  a 
Combination  of  Mach  1.4  and  Mach  2.0  Nozzle  Blocks 


At  Mach  1.5,  a  pressure  difference  of  approximately  7  percent  results 
in  a  Mach  number  change  of  only  0.05.  Therefore,  the  calculated  Mach  num- 
bers are  accurate  within  ±  0.05,  and  a  predicted  Mach  number  of  M^  =  1.50 
±  0.05  at  the  probe  tip  appears  reasonable.  Further  evidence  supporting  a 
value  of  Ma,  =  1.5  can  be  seen  in  Figure  4-3,  which  is  a  Schlieren  photo- 
graph of  the  model  in  the  wind  tunnel  test  section  during  one  of  the  test 
runs.  The  Mach  waves  emanating  from  the  top  and  bottom  of  the  test  section 
are  oriented  at  38  and  46  degrees,  respectively.  At  Mach  angle  y  =  46 
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degrees,  the  value  of  M*,  is  1.38.  The  probe  conical  nose  shock  angles  are 
<j>  =  38  (top)  and  A  =  46  (bottom),  indicating  that  the  probe  is  at  angle- 
of -attack.  A  rough  estimate  of  the  angle-of -attack  is  given  by: 


*2  "  *1 
a  =  — .  (4-1) 


With  <j>  =38  and  <J>  =46,  the  value  of  a  is  4  degrees.  An  accurate  value 
for  a  would  require  further  study.  An  angle-of-attack  between  4  and  6  de- 
grees was  predicted  theoretically.  Since  all  shocks  in  the  test  section 
are  curved  due  to  Mach  gradient,  all  Mach/shock  wave  angles  were  measured 
locally,  i.e.,  in  close  proximity  to  the  physical  object  generating  them. 
In  addition,  Figure  4-3  shows  a  shock  wave  just  forward  of  the  probe  tip; 
this  shock  wave  was  predicted  theoretically,  and  is  shown  in  Figure  3-4. 

Due  to  the  Mach  number  gradient  and  confused  flow  in  the  test  section, 
the  pressure  measured  by  the  probe  is  not  accurate,  and  hence  is  not  rele- 
vant. Table  4.5,  which  presents  sample  wind  tunnel  test  data  with  the 
asymmetric  nozzle  blocks,  is  included  for  completeness  only.  The  unortho- 
dox wind  tunnel  tests  did  prove,  however,  that  the  missile  nose  shock  will 
not  detach  at  Mach  numbers  above  1.5. 

Figure  4-4  is  a  graph  of  distance  versus  theoretical  Mach  number  and 
pressure  (PSIA)  in  the  wind  tunnel  test  section.  The  theoretical  pressure 
is  indicated  by  the  character  "0".  Pressures  measured  by  the  probe,  and 
the  Mach  numbers  calculated  from  the  Mach  waves  emanating  from  the  wind 
tunnel  walls,  are  also  indicated.  Moreover,  Mach  waves  at  several  points 
in  the  plane  of  the  probe  nose  were  measured  and  recorded  on  the  graph. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  following  conclusions  result  from  the  static  pressure  probe  study: 

(1)  The  static  pressure  probe  of  Ref.  1,  with  8  =  10  degrees  and 
oj  =  3  degrees,  will,  as  an  integral  part  of  the  missile  nose, 
perform  within  altitude  specifications  of  25,000  feet  ±  2,000 
feet  at  Mach  2.0  at  angles-of-attack  from  -  8  to  +  8  degrees. 

(2)  At  Mach  2.0,  within  an  angle-of -attack  range  from  0  to  6  degrees, 
the  probe  will  measure  free-stream  static  pressure  within  4  per- 
cent. A  4  percent  error  in  measurement  is  equivalent  to  an  alti- 
tude error  of  900  feet. 

(3)  The  missile  nose  shock  remains  downstream  of  the  probe  pressure 
ports  for  flight  Mach  numbers  above  1.5. 

(4)  Since  the  theory  predicted  probe  performance  within  2.6  percent 
at  Mach  2.0,  and  since  the  theory  predicted  satisfactory  per- 
formance at  Mach  1.4,  the  static  pressure  probe  should  perform 
within  specifications  at  Mach  1.4. 

B.  RECOMMENDATIONS 

The  minimum  operating  Mach  number  for  the  static  pressure  probe  de- 
pends on  the  position  of  the  missile  nose  shock  wave.  At  some  lower  Mach 
number,  the  missile  nose  shock  wave  moves  upstream  of  the  pressure  ports. 
M  is  defined  as  the  minimum  Mach  number  for  satisfactory  probe  operation. 
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The  potential  use  of  the  static  pressure  probe  as  an  altitude-measuring 

device  is  limited  bv  the  minimum  Mach  number,  M  ,  and  by  excursions  in  the 

m 

angle-of-attack.  Recommended  further  studies  are  as  follows: 

(1)  Replace  the  missing  supersonic  wind  tunnel  plenum  chamber  screen, 
design  and  manufacture  a  set  of  Mach  1.4  nozzle  blocks  with  a 
larger  second  throat,  and  test  the  probe.  Determine  whether 

M  <  1.4  or  M  >  1.4. 

m  m 

(2)  Once  the  actual  flight  Mach  number  has  been  determined,  recalcu- 
late the  theoretical  distribution,  if  required,  and  reposition 
the  probe  pressure  ports  to  remove  the  error  due  to  Mach  number 
range--an  error  presently  inherent  in  probe  design. 
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4.0  5.0  6.0 

Mach  number,  Mx 


7.o        ao 


9.0 


Figure  2-2     Shock  Angle  Versus  Mach  Number  for  Various  Cone  Angles 
(Reproduced  from  Kuethe  and  Chow   [2]) 
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Leeward  Side 


Windward  Side 


Figure  2-8     Static  Pressure  Probe  Cross  Section 
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APPENDIX  A:   PRESSURE  DISTRIBUTION  PROGRAM  LISTING  AND  SAMPLE  DATA  RUNS 

Note:      In  this  typed  listing,  and  in  the  sample  data  runs,  the  character  "°" 
signifies  "greater  than,"  and  "±"  signifies  "less  than."  The  ampersand  (§) 
denotes  exponentiation. 

00001  LPRINT'SUPERSONIC  FLOW  PRESSURE  COEFFICIENTS  ABOUT  A  YAWED/' 

00002  LPR I  NT 'SLENDER,  POINTED,  BODY  OF  REVOLUTION.  PROGRAM  BY" 

00003  LPRINT"J.R.SCHONBERGER  AND  K.D.TILLOTSON.  BASED  ON  " 

00004  LPRINT'CALCULATION  PROCEDURES  BY  PROF.  T.H.GAWAIN." 

00005  LPRINT'NAVAL  POSTGRADUATE  SCHOOL,  MONTEREY  CA.  FALL  1980" 

00006  CLEAR 

00010  CLS:PI=3. 14159 

00020  DIM  A(20,20),B(20,20),C(20,20),CP(20),D(20,20),FPC20),G(20),GA(2 
0),GBC20),GCC20),GNC20),GO(20),GPC20),IXC20),JOC20),JR(20),JX(20 
),RC20),RPC20),T(20,20),X(20),PRC20) 


00025  LPRINT'ENTER  XI"; 
00027  LPRINT"ENTER  X2"; 
00029  LPRINT'ENTER  X3"; 
00031  LPRINT'ENTER  Bl"; 
00033  LPRINT'ENTER  B2"; 


INPUT  XlrLPRINT  XI 

INPUT  X2ILPRINT  X2 

INPUT  X3:LPRINT  X3 

INPUT  BULPRINT  Bl 

INPUT  B2:LPRINT  B2 
00035  LPRINT'ENTER  COEFF  A"; : INPUT  AA:LPRINT  AA 
00037  LPRINT"ENTER  COEFF  B"; : INPUT  BBrLPRINT  BB 
00039  LPRINTMENTER  COEFF  C"; : INPUT  CC:LPRINT  CC 
00041  LPRINT'ENTER  RX2M; : INPUT  RXrLPRINT  RX 
00043  LPRINT'ENTER  FREE  STREAM  MACH  NBR";:INPUT  M:LPRINT  M 
00045  MU=ATNCl/M/SQRC-l/M::l/M+l)) 
00047  RY=RX+TAN(MU)"(X3-X2) 

00049  LPRINT"ENTER  NBR  OF  SUB INTERVALS" ; : INPUT  N:LPRINT  N:J=N-1 
00051  LPRINT'ENTER  LENGTH  OF  SUBINTERVAL"; : INPUT  H:LPRINT  H 
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00053    FOR   Z=1T0J 

00055  XCZ)=Z"H/C1-TAN(B1)/TANCMU)) 

00057  RCZ)=X(Z)«"TAN(B1) 

00059  RP(Z)=TAN(B1) 

00061  IF  XCZ)±=X1  THEN  NEXT  Z  ELSE  63 

00063  FOR  Y=ZTOJ 

00065  R=1-AA/CTANCMU))S2 

00067  S=-1"C2"Y"H+BB/CTANCMU))S2) 

00069  T=-1"CC/CTAN(MU))S2+(Y"H)&2 

00071  XCY)=C-S+SQRCSS2-V:R"T))/(2"R) 

00073  R(Y)=SQR(AA"X(Y)S2+BB"X(Y)+CC) 

00075  RPCY)=C2J'AA"X(Y)+BB)/C2-R(Y)) 

00077  IF  X(Y)±=X2  THEN  NEXT  Y  ELSE  79 

00079  FOR  Z=YTOJ 

00081  XCZ)=CCRX-X2"TANCB2))/TAN(MU)+Z"H)/(1-TAN(B2)/TANCMU)) 

00083  R(Z)=TAN(B2)"X(Z)+RX-TAN(B2)"X2 

00085  RP(Z)=TAN(B2) 

00087  IF  XCZ)±=X3  THEN  NEXT  Z  ELSE  89 

00089  FOR  Y=ZTOJ 

00091  X(Y)=RY/TAN(MU)+Y"H 

00093  RCY)=RY 

00095  RP(Y)=0 

00097  NEXT  Y 

00099  LPRINT  TABCl)"J";TABCl0)"XCJ)M;TABC20)"RCJ)";TAB(30)nRPCJ)n 

00101  FOR  Z=1T0J 

00103  LPRINT  TAB(1)Z;TABC10)XCZ);TABC20)RCZ);TABC30)RPCZ) 
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00110  IF  Z=10  THEN  INPUT"PRESS  ENTER  TO  CONTINUE";N$ 

00120  NEXT  Z 

00190  B=SQRCMS2-1) 

00192  ON  ERROR  GOTO  1200 

00195  U=0 

00200  FOR  Z=1T0J 

00205  U=U+1 

00210  FOR  Y=0TOU 

00220  TCZ,Y)=SQRCCXCZ)-Y"H)S2-CB"R(Z))&2) 

00230  NEXT  Y 

00231  NEXT  Z 
00235  U=0 
00240  FOR  Z=1T0J 
00245  U=U+1 
00250  FOR  Y=1T0U 

00260  ACZ,Y)=LOGCCCXCZ)-CY-l)"H)+T(Z,Y-l))/CCXCZ)-Y:'H)+TCZ,Y))) 

00270  NEXT  Y 

00275  NEXT  Z 

00280  FOR  Z=1T0J 

00290  FOR  Y=1T0J 

00300  B(Z,Y)=CT(Z,Y-1)-TCZ,Y))/RCZ) 

00310  NEXT  Y 

00311  NEXT  Z 
00320  FOR  Z=1T0J 
00330  FOR  Y=1T0J 

00340  C(Z,Y)=B(Z,Y)-RP(Z)"A(Z,Y) 
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00350  NEXT  Y 

00351  NEXT  Z 
00355  U=0 

00360  FPC1)=RPC1)/C(1,1) 
00370  FOR  Z=2  TO  J 
00380  U=U+1 
00390  FOR  Y=1T0U 
00400  W=W+C(Z,Y)::FP(Y) 

00410  NEXT  Y 

00420  FPCZ)=(RP(Z)-W)/C(Z,Z) 

00430  W=0 

00440  NEXT  Z 

00450  FOR  Z=1T0J 

00460  FOR  Y=1T0J 

00470  W=W+A(Z,Y>'FP(Y) 

00480  NEXT  Y 

00490  IX(Z)=W 

00500  W=0 

00510  NEXT  Z 

00520  FOR  Z=1T0J 

00540  NEXT 

00550  FOR  Z^ITOJ 

00560  GOCZ)=Cl-Cl+RPCZ)S2)"Cl-IXCZ))S2+CM"IX(Z))S2) 

00565  NEXT  Z 

00610  FOR  Z=l  TO  J 

00620  FOR  Y=l  TO  Z 
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00630  DCZ,Y)=(CXCZ)-CY-1)"H)"TCZ,Y-1)-CX(Z)-Y"H)"TCZ/Y))/C2"RCZ)S2) 

00640  NEXT  Y 

00650  NEXT  Z 

00660  FOR  Z=l  TO  J 

00670  FOR  Y=l  TO  Z 

00680  D(Z,Y)=-DCZ,Y)+BS2"ACZ,Y)+(XCZ)-(Y-1)"H)"BCZ/Y)/RCZ) 

00690  NEXT  Y 

00700  NEXT  Z 

00710  FOR  Z=1T0J 

00720  FOR  Y=1T0Z 

00730  DCZ,Y)=(D(Z,Y)+BS2J*A(Z,Y)+RP(Z):sB(Z,Y)) 

00740  NEXT  Y 

00750  NEXT  Z 

00760  GC0)=0:GPCD=1/DC1,1):U=0 

00770  FOR  Z=2TOJ 

00780  W=0:V=0 

00790  U=U+GP(Z-1) 

00800  G(Z-1)=H"U 

00810  FOR  Y=1T0Z 

00820  W=W+D(Z,Y-1)"GP(Y-1) 

00830  V=V+B(Z,Y)"GCY-1) 

00840  NEXT  Y 

00850  GP(Z)=C1-V/RCZ)-W)/DCZ,Z) 

00860  NEXT  Z 

00870  FOR  Z=1T0J 

00880  U=0:V=0 
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00890  FOR  Y=1T0Z 

00900  U=U+B(Z,Y)"GP(Y) 

00910  V=V+A(Z,Y)::GP(Y) 

00920  NEXT  Y 

00930  JX(Z)=U 

009^0  JRCZ)=VICBS2 

00950  JO(Z)=l-JR(Z)-RP(Z)"JX(Z) 

00970  NEXT  Z 

00990  FOR  Z=1T0J 

01010  GN(Z)=((l+RP(Z)S2)"(l-IX(Z))-fM&2:*IX(Z))5\JX(Z) 

01020  GA(Z)=-1"MS2"IX(Z)&2+(1+RP(Z)S2)"(1-IX(Z))&2 

01030  GB(Z)=(MS2-1-RP(Z)&2)"JX(Z)S2 

010<+0  GC(Z)=(l+JO(Z))&2 

01050  NEXT  Z 

01055  CLS 

01060  LPRINT'ENTER  THE  ANGLE  OF  ATTACK  A  (IN  RADIANS)"; : INPUT  A:LPRINT 
A 

01070  LPRINT"ENTER  THE  RADIAL  ANGLE  THETA  (IN  RADIANS)"; : INPUT  T:LPRIN 
T  T 

01071  LPRINT"ENTER  TEMP  (KELVIN)"; : INPUT  TElLPRINT  TE 

01072  LPRINT"ENTER  PRESSURE  (N/MS2)"; : INPUT  PE:LPRINT  PE 

01073  LPRINr'ENTER  DENSITY  (KG/MS3)"; : INPUT  POrLPRINT  PO 

01080  FOR  Z=1T0J 

01090  CP(Z)=GO(Z)-2"GN(Z)"SIN(A)"COS(A)"COS(T)+(GA(Z)4GB(Z)"(COS(T))S2 
-GC(Z)"(SIN(T))S2)"(SIN(A))S2 

01095  PR(Z)=l+200.9i+l"PO:cTE"MS2"CP(Z)/PE 

01100  NEXT  Z 
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01120  LPRINTTABCl)"Ju;TAB(10)"XCJ)";TAB(20)MCP(X)n;TAB(32)"PRESSURE  RA 
TIO" 

01130  FOR  Z=1T0J 

01140  LPRINTTAB(1)Z;TAB(10)X(Z);TAB(20)CP(Z);TAB(32)PR(Z) 

01160  NEXT  Z 

01170  LPRINTIF  YOU  DESIRE  TO  CHANGE  A  AND  THETA  ENTER  Y  ELSE  NM;:INPU 
T  N$:LPRINT  N$ 

01180  IF  N$="N"  THEN  END  ELSE  GOTO1055 
01200  X(Z)=X(Z)+X(Z)::.001 
01210  RESUME 
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Sample  Data  Run  for  Maah  2.0  at  Angle  of  Attack  a  =  0  and  10  Degrees 


SUPERSONIC  FLOW  PRESSURE  COEFFICIENTS  ABOUT  A  YAWED, 
SLENDER,  POINTED,  BODY  OF  REVOLUTION.  PROGRAM  BY 
J.R.SCHONBERGER  AND  K.D.TILLOTSON.  BASED  ON 
CALCULATION  PROCEDURES  BY  PROF.  T.H.GAWAIN. 
NAVAL  POSTGRADUATE  SCHOOL,  MONTEREY  CA.  FALL  1980 


ENTER 

XI 

100 

ENTER 

X2 

260 

ENTER 

X3 

400 

ENTER 

Bl 

.17453 

ENTER 

B2 

.0523599 

ENTER 

COEFF  A-9.26731E-03 

ENTER 

COEFF  B  8.36762 

ENTER 

COEFF  C-402.881 

ENTER 

RX2 

33.856 

ENTER 

FREE  STREAM  MACH  NBR  2 

ENTER 

NBR 

OF  SUB INTERVALS  20 

ENTER 

LENGTH  OF  SUB INTERVAL  15 

J 

XCJ) 

RCJ) 

RPCJ) 

1 

21.5952 

3.80776 

.176324 

2 

43.1905 

7.61551 

.176324 

3 

64.7857 

11.4233 

.176324 

4 

86.3809 

15.231 

.176324 

5 

109.828 

20.1081 

.157449 

6 

129.756 

22.9529 

.129889 

7 

148.696 

25.2279 

.111218 

8 

166.989 

27.1293 

.0971743 

9 

184.81 

28.7579 

.0859282 

10 

202.262 

30.1736 

.0765366 

11 

219.414 

31.4156 

.0684513 

12 

236.311 

32.511 

.0613283 

13 

252.989 

33.4798 

.0549372 

14 

269.503 

34.354 

.0524078 

15 

286 

35.2186 

.0524078 

16 

302.498 

36.0832 

.0524078 

17 

318.996 

36.9478 

.0524078 

18 

335.493 

37.8124 

.0524078 

19 

351.991 

38.677 

.0524078 

ENTER 

THE 

ANGLE  OF  , 

ATTACK  A  CIN  RADIANS)  0 

ENTER 

THE 

RADIAL  ANGLE  THETA 

CIN  RADIANS)  0 

ENTER 

TEMP  (KELVIN) 

288.15 

ENTER 

PRESSURE  (N/MS2)  101325 

ENTER 

DENSITY  (KG/M&3)  1.225 
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J  X(J)  CP(X) 

1  21.5952  .114208 

2  43.1905  .114208 

3  64.7857  .114208 

4  86.3809  .114208 

5  109.938  .0989367 

6  129.756  .0651711 

7  148.696  .0455148 

8  167.156  .0324173 

9  184.81  .0229932 

10  202.464  .015794 

11  219.414  .0100775 

12  236.311  5.30874E 

13  252.989  1.29649E 

14  269.503  1.17926E 

15  286  3.38801E 

16  302.498  5.17499E 

17  318.996  6.62418E 

18  335.493  7.80804E 

19  352.343  8.76572E 


IF  YOU  DESIRE  TO  CHANGE  A  AND 
ENTER  THE  ANGLE  OF  ATTACK  A  ( 
ENTER  THE  RADIAL  ANGLE  THETA 
ENTER  TEMP  (KELVIN)  288.15 
ENTER  PRESSURE  (N/MS2)  101325 
ENTER  DENSITY  (KG/MS 3)  1.225 
J        X(J)      CP(X) 

1  21.5952   .132234 

2  43.1905   .132234 

3  64.7857   .132234 

4  86.3809   .132234 

5  109.938   .117314 

6  129.756   .0821772 

7  148.696   .0616945 

8  167.156   .0482622 

9  184.81    .0377925 

10  202.464   .0303393 

11  219.414   .0236319 

12  236.311   .018247 

13  252.989   .0136563 

14  269.503   .0130074 

15  286       .0147623 

16  302.498   .0161278 

17  318.996   .0171903 

18  335.493   .0180226 

19  352.343   .0190157 
IF  YOU  DESIRE  TO  CHANGE  A  AND 


PRESSURE  RATIO 

1.31979 

1.31979 

1.31979 

1.31979 

1.27703 

1.18248 

1.12744 

1.09077 

1.06438 

1.04422 

1.02822 
-03  1.01486 
-03  1.00363 
-03  1.0033 
-03  1.00949 
-03  1.01449 
-03  1.01855 
-03  1.02186 
-03  1.02454 

THETA  ENTER  Y  ELSE  NY 
IN  RADIANS)  .0349 
(IN  RADIANS)  3.1416 


PRESSURE  RATIO 
1.37026 
1.37026 
1.37026 
1.37026 
1.32849 
1.2301 
1.17275 
1.13514 
1.10582 
1.08495 
1.06617 
1.05109 
1.03824 
1.03642 
1.04134 
1.04516 
1.04813 
1.05046 
1.05325 
THETA  ENTER  Y  ELSE  NN 
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APPENDIX  B:   COMPUTER  PROGRAM  AND  SAMPLE  DATA  RUN 
FOR  CALCULATING  PROBE  PLENUM  PRESSURE 


Program 

00010  CLS 

00020  INPUT  "ENTER  PI"; PI 

00030  INPUT  "ENTER  P2";P2 

00040  INPUT  "ENTER  P3";P3 

00050  PRINT  "T=";T/,PP=";PP,"CASE=";CA 

00060  INPUT  "ENTER  YOUR  GUESS  FOR  PP";PP 

00070  A=P1-PP 

00080  B=P2-PP 

00090  C=P3-PP 

00100  IF  B<=0  AND  C<=0  THEN  200  (Case  3) 

00110  IF  C<=0  THEN  300  (Case  2) 

00120  IF  B<=0  THEN  400  (Case  1) 

00130  PRINT"SOMETHING  WRONG"; :STOP 

00200  T=SQR(2"ABS(A))-2»SQR(2"ABS(B))-SQR(2!SABS(C)) 

00210  CA=3 

00220  GOTO  50 

00300  T=SQR(2"ABS(A))+2"SQR(2"ABS(B))-SQR(2"ABS(C)) 

00310  CA=2 

00320  GOTO  50 

00400  T=SQR(2"ABS(A))-2"SQRC2"ABS(B))+SQR(2"ABS(C)) 

00410  CA=1 
00420  GOTO  50 
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Sample  Data  Run   (Reproduced  from  CRT  display) 


ENTER  PI  71.019237 

ENTER  P2  7.96525 

ENTER  P3  7.964412 

T=  0    PP=  0    CASE=  0 

ENTER  YOUR  GUESS  FOR  PP  ?.977 

T=  .174618    PP=  .977    CASE=  3 

ENTER  YOUR  GUESS  FOR  PP  ?.970 

T=  .0131546    PP=  .970    CASE=  3 

ENTER  YOUR  GUESS  FOR  PP  7.971 

T=-. 0186603    PP=  .971    CASE=  3 

ENTER  YOUR  GUESS  FOR  PP  7.9705 

T=-3.07509E-03    PP=  .9705    CASE=  3 

ENTER  YOUR  GUESS  FOR  PP  7.9704 

T=  1.14933E-04    PP=  .9704    CASE=  3 

ENTER  YOUR  GUESS  FOR  PP  7.970404 

T=-1.11237E-05    PP=  .970404    CASE=  3 
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APPENDIX  C:   SEQUENTIAL  LISTING  OF  TESTING 

As  stated  in  Chapter  3,  Appendix  C  is  a  sequential  listing  of  testing 
and  corrective  action  taken  between  testing  attempts,  and  a  discussion  of 
the  three  most  likely  reasons  for  the  failure  of  the  wind  tunnel  to  start 
at  Mach  1.4.  Testing  attempts,  and  the  corrective  actions  taken  between 
testing  attempts,  were  as  follows: 

(1)  Test  Number  One  was  conducted  with  the  Mach  1.4  nozzle  blocks 
installed  and  the  model  mounted  in  the  wind  tunnel  test  section.  Super- 
sonic flow  could  not  be  established  in  the  test  section.  Test  section 
blockage  was  considered  to  be  the  most  probable  cause.  A  higher  percent- 
age, 28%,  of  test  section  blockage  is  allowed  at  Mach  2;  therefore,  the 
decision  was  made  to  complete  Mach  2  testing  prior  to  reducing  model  and 
mounting  system  frontal  area.  The  Mach  1.4  nozzle  blocks  were  removed. 
The  wind  tunnel  plenum  chamber  screen  was  discovered  to  be  torn,  and  was 
hanging  loose  in  the  plenum  chamber.  The  two  functions  of  the  screen  are: 
(a)  to  obtain  a  flat  velocity  profile  (flat  stagnation  pressure  profile) , 
and  (b)  to  generate  small-scale  turbulence  which  dissipates  prior  to  ar- 
rival at  the  test  section.  Sufficient  manpower  to  replace  the  plenum  cham- 
ber screen  was  not  available.  Therefore,  the  screen  was  removed. 

(2)  The  Mach  2.0  nozzle  blocks  were  installed.  To  insure  that  the 
flow  in  the  wind  tunnel  test  section  was  uniform,  an  additional  pressure 
tap  was  installed  in  the  wind  tunnel  wall  opposite  one  of  the  existing 
pressure  taps.  Testing  showed  no  significant  pressure  difference  (approxi- 
mately 0.5  percent)  between  the  two  pressure  taps.  Therefore,  the  test 
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section  velocity  profile  was  inferred  to  be  reasonably  uniform.  Twenty- 
three  Mach  2.0  wind  tunnel  tests  were  conducted  with  supersonic  flow 
achieved  in  the  test  section  (25  to  30  PSIG  plenum  chamber  pressure) . 

(3)  The  Mach  2.0  nozzle  blocks  were  removed,  and  the  Mach  1.4  noz- 
zle blocks  were  reinstalled.  The  model  was  mounted  in  the  wind  tunnel, 
and  testing  was  attempted.  Supersonic  flow  in  the  test  section  was  not 
achieved. 

(4)  The  model  frontal  area  was  reduced  by  14  percent  (yielding  5.4 
percent  test  section  blockage)  by  machining  0.050  inch  from  the  top  and 
bottom  of  the  missile  nose  ogive.  Supersonic  flow  still  could  not  be 
achieved. 

(5)  Model  frontal  area  then  was  reduced  by  an  additional  30  percent 
(yielding  3.5  percent  test  section  blockage)  by  machining  an  additional 
0.060  inch  from  the  top  and  bottom  of  the  missile  nose  ogive  and  0.015 
inch  from  the  top  and  bottom  of  the  double  edge  airfoil.  After  the  reduc- 
tion in  blockage,  supersonic  flow  could  not  be  achieved. 

(6)  Four  boundary- layer  bleed  valves  were  installed  in  both  wind 
tunnel  walls,  just  forward  of  the  nozzle  second  throat,  to  decrease  the 
required  mass  flow  rate  through  the  second  throat.  A  complete  description 
of  the  bleed  valve  specifications  and  function  can  be  found  in  Habel  [8] . 
Once  again,  supersonic  flow  could  not  be  achieved. 

(7)  The  bleed  valves  and  the  model  were  removed  from  the  wind  tun- 
nel. Attempts  to  obtain  supersonic  flow  in  the  wind  tunnel  test  section 
failed  with  no  wind  tunnel  blockage. 

(8)  A  fine-mesh  screen  was  installed  in  the  wind  tunnel  just  for- 
ward of  the  nozzle  first  throat.  With  the  screen,  but  without  the  model, 
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supersonic  flow  could  not  be  achieved.  Figure  A-l  is  a  photograph  of  the 
shock  that  was  observed  in  the  wind  tunnel  test  section  between  25  and  65 
PSIG  plenum  chamber  pressure. 

The  three  most  probable  causes  of  the  failure  to  obtain  supersonic 
flow  were:  (a)  that  the  model  was  too  large,  resulting  in  choked  flow;  (b) 
that  the  wind  tunnel  plenum  chamber  screen  was  missing,  causing  an  uneven 
velocity  profile;  and  (c)  that  the  Mach  1.4  nozzle  block  second  throat  was 
too  small,  resulting  in  choked  flow.  These  factors,  of  course,  can  operate 
in  combination  as  well  as  singly. 

Model  size  was  ruled  out  as  a  cause  of  failure  due  to  two  considera- 
tions. First,  after  reduction  of  the  model's  frontal  area,  the  actual  wind 
tunnel  blockage,  3.5  percent,  was  well  below  the  6.4  percent  theoretical 
allowed  wind  tunnel  blockage.  Second,  supersonic  flow  could  not  be  achieved 
even  without  the  model  in  the  wind  tunnel.  The  wind  tunnel  plenum  chamber 
screen  and  the  nozzle  block  second  throat,  however,  remained  as  probable 
causes  of  testing  failure. 

The  inside  dimensions  of  the  Mach  1.4  nozzle  blocks  were  measured, 
and  the  areas  of  the  first  and  second  throats  were  14.36  and  15.61  square 
inches,  respectively.  The  ratio  of  second  throat  area  to  first  throat  area 
was  thus  1.09.  Using  the  calculation  method  presented  in  Liepman  and  Roshko 
[6],  the  theoretical  minimum  area  ratio  allowable  is  1.04.  Therefore,  the 
Mach  1.4  nozzle  block  design  allows  a  safety  factor  of  approximately  five 
percent  for  the  boundary  layer  in  the  second  throat. 

Since  there  is  less  constriction  in  the  Mach  1.4  nozzle  than  in  the 
Mach  2.0  nozzle,  an  uneven  velocity  profile  would  affect  more  adversely 
the  flow  in  the  Mach  1.4  nozzle.  The  adverse  effect  of  the  uneven  velocity 
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profile  may  have  enlarged  the  boundary  layer  in  the  second  throat,  thus 
preventing  supersonic  flow. 
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